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' Abstract. High resolution echelle spectra of PW And (HD 1405) have been taken during eight observing runs 

from 1999 to 2002. The detailed analysis of the spectra allow us to determine its spectral type (K2V), mean 
heliocentric radial velocity (Vhei = —11.15 km s~^) rotational velocity (usini = 22.6 km s~^), and equivalent 
> I width of the lithium line A6707.8 A (£;W(Li I) = 273 mA). The kinematic (Galactic Velocity {U,V,W)) confirms 

its membership of the Local Association moving group, in agreement with the age (30 to 80 Myrs) inferred 
from the color magnitude diagram and the lithium equivalent width. Photospheric activity (presence of cool 
spots that disturb the profiles of the photospheric lines) has been detected as changes in the the bisectors of the 
cross correlation function (CCF) resulting of cross-correlate the spectra of PW And with the spectrum of a non 
' active star of similar spectral type. These variations of the CCF bisectors are related to the variations in the 

, measured radial velocities and are modulated with a period similar to the photometric period of the star. At the 

same time, chromospheric activity has been analyzed, using the spectral subtraction technique and simultaneous 
spectroscopic observations of the Ha, H/3, Na I Di and D2, He I D3, Mg I b triplet, Ca 11 H&K, and Ca 11 infrared 
triplet lines. A flare was observed during the last observing run of 2001, showing an enhancement in the observed 



(N 

o 
o 



Q ' chromospheric lines. A less powerful flare was observed on 2002 August 23. The variations of the chromospheric 

activity indicators seem to be related to the photospheric activity. A correlation between radial velocity, changes 
C/3 ■ in the CCF bisectors and equivalent width of different chromospheric lines is observed with a different behaviour 



between epochs 1999, 2001 and 2002. 
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with emission. No clear luminosity classification has been 
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proposed but Hooten & Hall (1990) suggest it is a main 
sequence star due to its short-period photometric varia- 
tions (P ^ 1.75 days). Christian ct al. (2001) identify the 
EUVE source EUVE J0018+308 as PW And. This star 
was classified as a neighborhood Pleiades-age K2 dwarf 
by Ambnistcr ct al. (1998). In this work wc conclude from 
its position in the color magnitude diagram that PW And 
could be a pre-main sequence star. This is consistent with 
the equivalent width of the lithium line A6707.8 A mea- 
sured by Wichmann et al. (2003) EW{U i) = 298 ml 
and the abundance of lithium obtained by Ambruster et 
al. (1994), logA/'(Li i) = 3.00 dcx. 

With a rotational velocity {v sin i) of 23.4 km mea- 
sured by Fekel (1997) this star is a good candidate for 
emission in chromosphcric lines like Balmcr and Ca n H 
& K lines. Bidelman (1985) confirms the existence of mod- 
erate emission in the Ca II H & K lines and Strassmeier 
ct al. (1988) include PW And as a candidate in the 
first catalogue of Chromospherically Active Binary Stars. 
Nevertheless, studies by Fehrenbach & Burnage (1982) 
and Griffin (1992) fail to confirm a binary nature as the 
radial velocity (14) remains constant with a value of —11.5 
and —10.5 km s~^ respectively. 

Space velocity components in addition to its luminos- 
ity class, rotational velocity, radial velocity as well as its 
chromosphcric activity all indicate that PW And is a can- 
didate for belonging to the Local Association stellar kine- 
matic group (for a detailed reference see Montes et al. 
2001a, b, hereafter MOla and MOlb). 

In this paper we are interested primarily in measur- 
ing the variations caused by magnetic activity in both the 
photospheric and chromosphcric lines in the optical spec- 
tral region. At the same time wc carry out a study of the 
fundamental properties of the star. This is accomplished 
using a large set of medium and high resolution echelle 
spectra that were obtained during eight observing runs 
from 1999 to 2002. This interval samples the rotational 
period during different epochs. The spectra analyzed here 
include all the optical chromosphcric activity indicators 
from the Ca ii H & K to Ca II IRT hues as well as the Li i 
A6707.8 line. 

In Sect. 2 we give the details of our observations 
and data reduction. Stellar parameters determined by us 
(spectral type, luminosity class, rotational and radial ve- 
locities and age estimation) are given in Sect. 3. In Sect. 4 
we present the study of the chromosphcric activity indica- 
tors for each observing run, while flare events are studied 
in more detail in Sect. 5. Photospheric variations are dis- 
cussed in Sect. 6, and Sect. 7 discusses how the photomet- 
ric variations are related to the chromosphcric variations. 
Finally in Sect. 8 we state the conclusions. 

2. Observations and data reduction 

The echelle spectra of P W And analyzed in this paper were 
obtained during eight observing runs which are described 
in detail in the below paragraphs. 



1) 2.2m-FOCES 1999/07 

This run took place on 24-29 July 1999 using the 
2.2 m telescope at the German Spanish Astronomical 
Observatory (CAHA) (Almerfa, Spain). The Fiber Optics 
Cassegrain Echelle Spectrograph (FOCES) (Pfeiffer et al. 
1998) was used with a 2048^ 15/i LORAL#lli CCD de- 
tector. The wavelength range covers from 3910 to 9075 A 
in 84 orders. The reciprocal dispersion ranges from 0.03 
to 0.07 A/pixel and the spectral resolution, measured as 
the full width at half maximum (FWHM) of the arc com- 
parison lines, ranges from 0.09 to 0.26 A. 

2) NOT-SOFIN 1999/11 

Observations taken on 26-27 November 1999 used the 
2.56 m Nordic Optical Telescope (NOT) located at the 
Observatorio del Roque de Los Muchachos (La Palma, 
Spain). The Soviet Finnish High Resolution Echelle 
Spectrograph (SOFIN) was used with an echelle grating 
(79 grooves/mm), Astromcd-3200 camera and a 1152x770 
pixel EEV P88200 CCD detector. The wavelength range 
covered is from 3525 to 10425 A in 44 orders. The recip- 
rocal dispersion ranges from 0.06 to 0.17 A/pixel and the 
spectral resolution (FWHM) from 0.14 to 0.32 A. 

3) INT-MUSICOS 2000/08 

Observations made on 5-11 August 2000 with the 2.5 m 
Isaac Newton Telescope (INT) at the Observatorio del 
Roque de Los Muchachos (La Palma, Spain) used the 
ESA-MUSICOS spectrograph. This is a fibre-fed cross- 
dispersed echelle spectrograph, built as a replica of 
the first MUSICOS spectrograph (Baudrand & Bohm 
1992) and developed as part of MUlti-SIte Continuous 
Spectroscopy (MUSICOS^) project. During this observing 
run, a 1024^ 24/z TEK5 CCD detector was used, yielding 
a wavelength coverage from 4430 A to 10225 A in 73 or- 
ders. The reciprocal dispersion ranges from 0.07 to 0.15 A 
and the spectral resolution (FWHM) from 0.16 to 0.30 A. 

4) NOT-SOFIN 2000/11 

This run took place on 10-13 November 2000 with the 
2.56 m Nordic Optical Telescope (NOT) located at the 
Observatorio del Roque de Los Muchachos (La Palma, 
Spain). The Soviet Finnish High Resolution Echelle 
Spectrograph (SOFIN) was again used in the same con- 
figuration as in run (2). 

5) 2.2m-FOCES 2001/09 

These observations were made on 21-24 September using 
the 2.2 m telescope at the German Spanish Astronomical 
Observatory (CAHA) (Almerfa, Spain). The Fibre Optics 
Cassegrain Echelle Spectrograph (FOCES) (Pfeiffer et al. 
1998) was used with a 2048^ 24/i Site#ld CCD detec- 
tor. The wavelength range covers from 3510 to 10700 A 
in 112 orders. The reciprocal dispersion ranges from 0.04 
to 0.13 A/pixel and the spectral resolution, determined 
as the full width at half maximum (FWHM) of the arc 
comparison lines, ranges from 0.08 to 0.35 A. 

6) SARG-TNG 2001/10 

Observations were taken on 10-11 October 2001 with 
the 3.5 m Telescopic Nazionale Galileo (TNG) located 

^ http://www.ucm.es/info/Astrof/MUSICOS.html 
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at the Observatorio del Roque de Los Muchachos (La 
Palma, Spain). The Spectrografo di Alta Resoluzione 
Gahleo (SARG) was used with an echelle grating 
(31.6 grooves/mm, R4) and the red cross-disperser (200 
grooves/mm), and a mosaic of two 2048X4096 CCDs with 
a 13.5/i pixel size. The wavelength range covers from 4960 
to 10110 A in 62 orders. The reciprocal dispersion ranges 
from 0.02 to 0.04 A/pixel and the spectral resolution 
(FWHM) from 0.08 to 0.17 A. 

7) HET-HRS 2001/12 - 2002/02 

Observations made on 19 December 2001 - 28 February 
2002 using the 9.2 m Hobby-Eberly Telescope at 
McDonald Observatory in Texas (USA) with the High 
Resolution Spectrograph (HRS). The detector is a mosaic 
of two Marconi 4096 x 2048 CCDs with 15 micron pix- 
els (Tull 1998). The wavelength range covers from 5040 to 
8775 A in 52 orders. The reciprocal dispersion ranges from 
0.06 to 0.11 A/pixel and the spectral resolution (FWHM) 
from 0.15 to 0.28 A. 

8) NOT-SOFIN 2002/08 

This took place on 21-29 August 2002 with the 
2.56 m Nordic Optical Telescope (NOT) located at the 
Observatorio del Roque de Los Muchachos (La Palma, 
Spain). The Soviet Finnish High Resolution Echelle 
Spectrograph (SOFIN) was used with an echelle grat- 
ing (79 grooves/mm), camera 2 and a 2052x2052 pixel 
PISKONOV Loral CCD detector. The wavelength range 
covers from 3525 to 10200 A in 42 orders. The recipro- 
cal dispersion ranges from 0.015 to 0.045 A/pixel and the 
spectral resolution (FWHM) from 0.05 to 0.15 A. 

The velocity resolution derived from the FWHM spec- 
tral resolution in the different observing runs ranges from 
4 to 10 km s"^ 

Non-active stars used as reference stars in the spectral 
subtraction technique and the radial velocity standards 
used in the radial velocity determinations are listed in 
Tables n] and 121 In Table we give the observing log. For 
each observation we list date, UT and the signal to noise 
ratio (S/N) obtained in the Ca ii H & K and Ha line 
regions. The stars marked with * in the table are standard 
radial velocity stars taken from Beavers et al. (1979). 

The spectra have been extracted using the standard 
reduction procedures in the IRAF^ package (bias subtrac- 
tion, flat-field division and optimal extraction of the spec- 
tra) . The wavelength calibrations were obtained by taking 
spectra of a Th-Ar lamp. This lamp gives a very rich emis- 
sion line spectrum over a wide wavelength range and is 
thus very useful for high resolution echelle spectra, includ- 
ing enough calibration lines in each order. The calibration 
was performed by fitting all the orders simultaneously (in 
both dispersion and cross dispersion directions) and the 
residuals obtained in the fit are typically lower than 0.04 
A. Finally, the spectra were normalized by using a cu- 



Table 1. Observing log 



^ IRAF is distributed by the National Optical Observatory, 
which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation. 



Name 


Day/Month 


UT 


S/N 


S/N 








(H&K) 


(Ua) 


2.2m-FOCES 1999 


PW And 


25 Jul 


04:09 


23 


123 




26 " 


00:57 


19 


110 




27 " 


02:25 


13 


85 




28 


00:57 


15 


76 




29 " 


03:01 


21 


122 




30 " 


01:57 


20 


130 


GJ 706 


25 Jul 


00:44 


37 


263 




25 " 


20:07 


35 


176 


/3 Oph * 


24 Jul 


23:48 


31 


182 




25 " 


19:51 


25 


222 




25 " 


19:55 


25 


268 




27 " 


20:02 


22 


164 




28 " 


19:57 


39 


332 




28 " 


20:01 


42 


327 




29 " 


19:46 


32 


314 


NOT-SOFIN 1999 


PW And 


26 Oct 


20:24 


25 


107 




27 " 


21:29 


20 


155 


HR 222 * 


26 Oct 


23:06 


27 


119 


INT-MUSICOS 2000 


PW And 


11 Aug 


03:36 




49 




14 " 


03:31 




69 


HR 222 * 


11 Aug 


03:20 




119 


NOT-SOFIN 2000 


PW And 


06 Oct 


01:35 


7 


129 




08 " 


00:53 


16 


89 




09 " 


00:03 


8 


101 


HR 222 * 


07 Oct 


01:44 


11 


234 


2.2m-FOCES 2001 


PW And 


23 Sop 


23:59 


14 


91 




25 " 


02:03 


16 


99 


GJ 706 


23 Sop 


19:52 


23 


149 


HR, 166 * 


23 Sop 


23:50 


30 


174 


TNG-SARG 2001 


PW And 


12 Oct 


04:05 




65 


HR 222 * 


12 Oct 


04:30 




136 


HET-HRS 2001-2002 


PW And 


20 Dec 


04:12 




239 




21 " 


03:37 




338 




22 " 


03:44 




252 




23 " 


04:00 




352 




26 " 


03:09 




184 




27 " 


03:07 




206 




28 '• 


03:31 




278 




29 " 


02:52 




67 




30 " 


02:59 




101 


HD 92588 * 


14 Feb 


09:09 




246 


NOT-SOFIN 2002 


PW And 


22 Aug 


04:16 


25 


173 




23 " 


04:24 


21 


161 




24 " 


04:39 


22 


145 




25 " 


01:55 


14 


105 




26 » 


05:15 


14 


108 




26 " 


05:37 


12 


77 




27 " 


05:06 


14 


106 




28 " 


04:19 


19 


151 




29 " 


02:37 


20 


138 


HR 222 * 


22 Aug 


06:05 


24 


273 




23 " 


06:06 


20 


302 




24 


06:13 


18 


204 




25 " 


02:31 


32 


277 




26 " 


05:53 


22 


218 




27 


06:04 


22 


219 




28 " 


05:19 


15 


158 




29 " 


03:04 


25 


182 
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Table 2. Stellar parameters 



Name 


HD 


SpT 


V-R 


B-V 


V sin i 
(km s-^) 


-Pphot 

(days) 


Vhol ± TV, 1 

(km s""^) 


U ±av 
(km s-i) 


y ± o-v 
(km s-i) 


W ±uvr 
(km s-i) 


PW And 


HD 1405 


K2V'' 


0.74*" 


1.04"^ 


22.6±0.4=' 


1.75'^ 


-11.15±0.05'' 


-5.42±0.33'' 


-28.69i0.63'' 


-17.94i0.74'' 


Ref. Stars 


HR 4182 * 
GJ 706 
HR 222 * 


HD 92588 
HD 166620 
HD 4628 


KlIV 
K2V 
K2V 



















" This work. 

^ Tlic V-R color index are obtained from the relation with spectral type given by Landolt-Bornstein (Schmidt-Kaler 1982). 
From Tycho-2 Catalogue (H0g et al. 2000). 

^ Photometric period calculated by Hooten & Hall (1990). 



bic spline polynomial fit to the observed continuum. The 
points used to fit the continuum in each order are chosen 
in regions judged free of spectral lines to reduce as much 
as possible the effects of line blending. 

3. Stellar parameters 

3.1. Spectral type and luminosity class 

To obtain an independent estimate of the spectral type 
of this star, we have compared our high resolution echelle 
spectra with that of inactive reference stars taken during 
the same observing run. The spectral type of the reference 
stars range from G to M, and the luminosity class from III 
to V. The analysis makes use of the program STARMOD de- 
veloped at Penn State University (Barden 1985) and mod- 
ified more recently by us. With this program a synthetic 
stellar spectrum is constructed from the artificially rota- 
tionally broadened, and radial-velocity shifted spectrum 
of an appropriate reference star. We obtained the best fit 
between observed and synthetic spectra when we use a 
K2-dwarf spectral type standard star. The uncertainty in 
this classification is of one spectral subtype as is typical 
in the MK spectral classification. In addition, in order to 
constrain the spectral type and luminosity class, several 
temperature-and gravity-sensitive photospheric lines have 
been studied in more detail and lead to a very similar 
result. 

3.2. Radial velocity and space motion 

Heliocentric radial velocities were determined using the 
cross-correlation technique. The spectra of the star were 
cross-correlated order by order, using the routine FXCOR 
in IRAF, against spectra of radial velocity standards of 
similar spectral types (the stars marked with * in Tables 
n]and[2Jl taken from Beavers et al. (1979). For each order, 
the velocity is derived from the position of the peak of the 
cross-correlation function (CCF) by fitting a Gaussian to 
the top of the function. Radial velocity errors are com- 
puted by FXCOR based on the fitted peak height and 
the antisymmetric noise as described by Tonry & Davis 
(1979). The radial velocities calculated for each order are 
weighted by their errors, and a mean value is obtained for 



each observation (see the values given in Table |7| for three 
representative observing runs). Orders including chromo- 
spheric features and prominent telluric lines were excluded 
when determining this mean velocity. Finally, a weighted 
mean radial velocity is determined using all the observa- 
tions from all the observing runs. In Table |2] we list this 
average heliocentric radial velocity (Vhei) and its associ- 
ated error (crVhoi)- 

As will be discussed in detail in Section 6 the pres- 
ence of star-spots in the photosphere of this active star 
can disturb the profile of the photospheric lines and in- 
duce variations in the peak of the CCF which cause vari- 
ations in the measured radial velocity. Differences up to 
3 km s^^ are measured in the spectra from the same ob- 
serving run (see Table [71). As this is larger than the sin- 
gle measurement error these variations are attributed to 
asymmetry variations of the absorption lines caused by 
the stellar activity. In addition, variations up to 6 km 
have been found in the radial velocities measured at dif- 
ferent epochs (from 1999 to 2002). These variations can 
be attributed to a) the effects on the line profiles of chang- 
ing levels of stellar activity (such as changes in the spot 
coverage) at different epochs, b) the use of different radial 
velocity standard which have different errors, or c) the 
changes in the zero-point offset of the different observing 
runs and spectrographs. However, the final mean radial 
velocity value we have obtained {vr = —11.15 km s""'^) 
is very similar to those given by Fehrenbach & Burnage 
(1982) (wr = -11.5 km s"^) and Griffin (1992) {v, = -10.4 
km s^^). Therefore, the radial velocity we report in this 
paper not only confirms the previously-measured radial 
velocities, but allow us to conclusively state that this 
star is not a chromospherically active binary, as is sug- 
gested by Strassmeier et al. (1988), since the component 
radial velocity amplitude of these kind of binaries are typ- 
ically much lager, ranging from 20 to 100 km s~^ (see 
Strassmeier et al. 1993). 

We have used this mean radial velocity together with 
the spectroscopic parallax (My — 6.4 from Landolt- 
Borstein (Schmidth-Kaler 1982) for a K2 dwarf gives a 
distance d — 27.54 pc) and proper motions from the 
Tycho-2 (H0g et al. 2000) Catalogue, to calculate the 
Galactic space-velocity components {U, V, W) in a right- 
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handed coordinated system (positive in the directions of 
the Galactic center, Galactic rotation, and the North 
Galactic Pole, respectively), as determined by MOla. The 
resultant new values are given in Table |21 Its position in 
the Galactic-velocity diagram confirms its membership in 
the Local Association moving group (see also MOlb and 
Lopez-Santiago et al. 2003). PW And is situated inside the 
subgroup B4 associated with the Pleiades cluster identi- 
fied by Asiain et al. (1999). 

3.3. Rotational velocity 

Fekel (1997) measured the rotational velocity (wsini) of 
PW And obtaining a value of 23. 4± 1.0 km s""'^, assuming 
spectral type G5. The mean value from the five Coravel 
traces is 21 ± 1 km s^^ and Griffin (1992) determined 
vsini = 21.5 km s~^. To determine an accurate rota- 
tional velocity of PW And, we again make use of the 
cross-correlation technique in our high resolution echelle 
spectra by using the routine FXCOR in IRAF. In each one 
of the observing runs the observed spectra of PW And 
were cross-correlated against the spectrum of the tem- 
plate star and the width (FWHM) of CCF determined. 
The calibration of this width to yield an estimation of 
V sin i is done by cross-correlating artificially broadened 
spectra of the template star with the original template 
star spectrum. Broadened template spectra are created 
for V sin i spanning the expected range of values by convo- 
lution with a theoretical rotational profile (Gray 1992) us- 
ing the program STARMOD. The resultant relationship be- 
tween V sin i and FWHM of the CCF was then fitted with 
a fourth-order polynomial. We have tested this method 
with stars of known rotational velocity and have obtained 
good agreement. The uncertainties on the vsini values 
obtained by this method have been calculated using the 
parameter R defined by Tonry & Davis (1979) as the ratio 
of the CCF height to the rms antisymmetric component. 
This parameter is computed by the IRAF task FXCOR and 
provides a measure of the signal to noise ratio of the CCF. 
Tonry & Davis (1979) show that errors in the FWHM of 
the CCF are proportional to {1 + R)~^ and Hartmann et 
al. (1986) and Rhode et al. (2001) found that the quan- 
tity ±i;sin j(l-|-i?)^^ provides a good estimate for the 90% 
confidence level of a vsini measurement. Thus, we have 
adopted ±wsini(l + R)~^ as a reasonable estimate of the 
uncertainties on our wsinz measurements. 

We have determined v sin i by this method in all the 
spectra available of PW And. Very similar vsini values 
are obtained in spectra on different nights and at different 
epochs. The vsini derived by this method is not affected 
by the asymmetry of the CCF caused by star-spots (see 
Section 6), since those changes are mainly in the peak 
of the CCF but the vsini determination is based on the 
FWHM of the CCF. The resulting error weighted mean 
for all the observing runs is 22.6 ± 0.4 km s~^, which is 
the value given in Table |21 Our value is closer to the value 
given by Fekel (1997). The minimum radius (-Rsini) cal- 
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1 

(B-V) 
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Fig. 1. Mv-{B — V) diagram. We have used pre-main se- 
quence isochrones from Siess et al. (2000). Continuous line 
corresponds to ZAMS and dotted lines are 3, 10, 30 and 
80 Myrs isochrones (from top to bottom). 

culated using our value of v sin i and the Pphot by Hooten 
& Hall (1990) is 0.78±0.11 R©, which is consistent with a 
K2 V spectral type. Adopting a radius R = 0.80 ±0.05 R© 
for a K2 dwarf from the Schmidt-Kaler (1982) tables, we 
obtain the inclination of the rotation axis i = 77 ± 9° for 
PW And. This is only an estimate, since the radius could 
be slightly different due to the possible pre-main sequence 
nature of this star inferred in the next subsection. 

3.4. Estimation of age: the color-magnitude diagram 

We estimate the age of PW And from the color- magnitude 
diagram {M^ vs. (B-V)). The value of has been 
adopted from Landolt-Bornstein tables (Schmidt-Kaler 
1982) for a dK2 type star. The color index B-V (listed 
in Table I^J has been obtained from Tycho-2 Catalogue 
(H0g et al. 2000) using transformations from Bt and Vt 
to Johnson indexes (section 1.3 from Hipparcos Catalogue, 
ESA 1997). Pre-main sequence isochrones from Siess et al. 
(2000) have been adopted in Fig. (see Lopez-Santiago et 
al. 2003). The continuous line corresponds to the ZAMS 
and the dotted lines are 3, 10, 30 and 80 Myrs isochrones 
respectively (from top to bottom) . In order to estimate the 
error bars for we have assumed a deviation of one spec- 
tral type from the dK2 suggested here (Kl to K3). Error 
in B~V index has been calculated taking into account er- 
rors in i?T and Vr given by Tycho-2 and their relation 
with B-V. PW And's position in the diagram over the 30 
Myrs isochrone (see Fig. ^ suggest that it is a pre-main 
sequence star. Even if we take into account the errors bars 
(or a large estimation of these) the star continue slightly 
above the ZAMS. 
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6700 6720 
Wavelength (A) 



Fig. 2. Spectrum of PW And in the region of Li i A6707.8 
line.The non-active dK2 star GJ 706 (HD 166620) is plot- 
ted as reference. 



o 
o 



Pleiades (Neuhauser et al. 1997) 

(Soderblom et al. 1993) 

- - IC 2602 (Montes et al. 2001b) 

- - - Hyades (Soderblom et al. 1993} 




FO 



GO 



KO 



MO 



Spectral Type 

Fig. 3. The figure shows lithium equivalent width 
{EW (Li i)) as a function of spectral type for members of 
the young open clusters: IC 2602 (dashed line), Pleiades 
(solid lines) and Hyades (dotted line) . Only the upper en- 
velope is plotted for IC 2602 and Hyades, while upper and 
lower envelopes are presented for Pleiades cluster. The star 
corresponds to the position of PW And in the diagram. 



3.5. Estimation of age: the Li i X6707.8 line 

In order to further constrain the age for PW And, the 
equivalent width of the resonance doublet of Li i at 
A6707.8 A has been measured. This line is included in our 
echelle spectra for all the observing runs and is an impor- 
tant diagnostic of age in late- type stars since the lithium is 
destroyed easily by thermonuclear reactions in the stellar 
interior. A representative spectrum in the Li i line region 
of the star is plotted in Fig. [5] As can be seen in this fig- 
ure the Li I absorption line is very prominent in PW And 
in comparison with a reference star of the same spectral 
type. At this spectral resolution and with the rotational 



velocity of the star the Li i line is blended with the nearby 
Fe I A6707.41 A line. We have corrected the total measured 
equivalent width, EW(Li i+Fe i), by subtracting the EW 
of Fe I calculated from the empirical relationship with {B- 
V) given by Soderblom et al. (1990). The error weighted 
mean value of the individual values of EW {Li i) mea- 
sured on different nights and over different observing runs 
is 266.6 ±0.4 mA. This value is a bit lower than that mea- 
sured by Wichmann et al. (2003) of 273 mA and corrected 
for the contribution of the Fe i A6707.41 A line. Our cor- 
rected iJW(Li l) value has been converted into abundances 
by means of the curves of growth computed by Pallavicini 
et al. (1987). Adopting an effective temperature (Tcff) of 
4900 K for PW And, we have obtained logiV(Li i) 3.4 
dex (on a scale where logiV(H) = 12.0). This estimation 
is considerably higher than the obtained by Ambruster et 
al. (1994) of log A^(Li i) = 3.0 dex, but these authors do 
not give any information about the temperature and the 
model they used. 

We compare the EW (Li i) of PW And with those of 
well-known young open clusters of different ages. In the 
EW(Li i) versus spectral type diagram (Fig. ^ we have 
over-plotted the upper envelope of the EW(Li i) of IC 
2602 (10-35 Myr), the Pleiades (78-125 Myr), and the 
Hyades (600 Myr), open clusters. For the Pleiades we 
adopt the upper envelope determined by Neuhauser et 
al. (1997) with data from Soderblom et al. (1993) and 
Garcia Lopez et al. (1994) and the lower envelope given 
by Soderblom et al. (1993). In the case of IC 2602 we have 
adopted the upper envelope determined in ALOlb. Finally 
for the Hyades (600 Myr) we have used the upper envelope 
adopted by Soderblom et al. (1993). 

The position of PW And close to the upper envelope 
of the Pleiades and under the upper envelope of IC 2602 
(see Fig. suggests an age similar to the stars in these 
clusters. Due to the large dispersion in EW{Li i) found 
for stars in young clusters, as in the case of IC 2602 (see 
MOlb), we estimate an age for PW And between 30 to 80 
Myrs, corresponding with the ages estimated for the IC 
2602 and the Pleiades clusters. This is consistent with the 
space motion obtained in section 3.2 and the age obtained 
by its position in the color-magnitude diagram (Fig. ^ . 

4. Chromospheric activity indicators 

The activity of PW And was previously studied by 
Bidelman (1985) who reported moderate Ca ii H & K 
emission and the Ha line in emission. Chromospheric and 
transition region UV emission fluxes have been reported 
by Ambruster et al. (1998) and Wood et al. (2000). In 
addition, it has been detected by the ROSAT-satellite as 
the 2RE J001820-f 305 source (Pye et al. 1995; Kreysing 
et al. 1995; Thomas et al. 1998; Rutledge et al. 2000), and 
by the EUVE-satellite as the EUVE JOO 18+309 source 
(Malina et al. 1994; Christian et al. 2001). 

The echelle spectra analyzed in this paper allow us to 
study the behaviour of different optical chromospheric ac- 
tivity indicators from the Ca ii H & K to the Ca ii IRT 
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Table 3. EW of the different chromospheric activity indicators 



EW(A) in the subtracted speetrum 



MJD 


Ca 11 














Ca 11 IRT 






K 


H 


He 


H(5 


H7 


H/3 


Ha 


A8498 


A8542 


A8662 


2.2m-FOCES 


1999 




















51384.1732 


1.78±0.09 


1.17±0.04 


0.31±0.01 


0.27±0.01 


0.29±0.01 


0.47±0.01 


1.39±0.04 


0.49±0.01 


0.63±0.01 


0.52±0.01 


51385.0401 


1.36±0.16 


0.99±0.06 


0.33±0.02 


0.28±0.01 


0.28±0.01 


0.45±0.01 


1.26±0.03 


0.51±0.01 


0.63±0.01 


0.56±0.01 


51386.1011 


1.98±0.27 


0.93±0.05 


0.32±0.02 


0.26±0.01 


0.29±0.01 


0.51±0.02 


1.35±0.03 


0.50±0.01 


0.66±0.01 


0.53±0.01 


51387.0396 








0.28±0.02 


0.27±0.01 


0.44±0.01 


1.23±0.03 


0.47±0.01 


0.62±0.01 


0.55±0.01 


51388.1258 


1.97±0.12 


0.82±0.02 


0.26±0.01 


0.24±0.01 


0.26±0.01 


0.48±0.01 


1.20±0.03 


0.47±0.01 


0.61±0.01 


0.52±0.01 


51389.0818 


2.18±0.32 


1.68±0.03 


0.33±0.04 






0.47±0.01 


1.36±0.04 


0.49±0.01 


0.68±0.01 


0.53±0.01 


NOT-SOFIN 


1999 




















51508.8690 


2.43±0.21 










0.71±0.02 


1.58±0.04 


0.58±0.01 


0.78±0.06 




51509.9022 


2.36±0.14 










0.71±0.04 


1.54±0.04 


0.53±0.01 


0.91±0.02 




INT-MUSICOS 2000 


51767.6572 










0.17±0.01 


0.49±0.01 


1.18±0.02 


0.49±0.01 


0.72±0.02 


0.59±0.02 


51770.6541 














0.84±0.05 


0.79±0.11 


0.66±0.06 


0.89±0.10 


NOT-SOFIN 


2000 




















51854.5731 


0.98±0.19 






0.20±0.04 


0.28±0.06 


0.51±0.02 


1.39±0.03 


0.49±0.01 


0.73±0.02 




51855.5441 


1.20±0.17 






0.25±0.03 


0.23±0.03 


0.54±0.02 


1.33±0.03 


0.47±0.01 


0.76±0.04 




51856.5440 


1.04±0.12 










0.60±0.02 


1.20±0.03 


0.45±0.01 


0.65±0.03 




51857.5090 


1.19±0.19 








0.30±0.05 


0.54±0.02 


1.48±0.04 


0.48±0.01 


0.67±0.02 




2.2m-FOCES 


2001 




















52176.4998 


1.79±0.07 


1.11±0.09 


0.39±0.03 


0.29±0.01 


0.26±0.01 


0.53±0.01 


1.25±0.03 


0.47±0.01 


0.70±0.01 


0.55±0.01 


52177.5860 


1.49±0.04 


1.02±0.06 


0.33±0.02 


0.25±0.01 


0.23±0.01 


0.47±0.01 


1.27±0.02 


0.43±0.01 


0.63±0.01 


0.51±0.01 



HET-HRS 2001-2002 



52263.6771 






1.26±0.02 


0.43±0.01 


0.71±0.02 


0.64±0.02 


52264.6541 






1.16±0.02 


0.50±0.02 


0.64±0.03 


0.55±0.02 


52265.6573 






1.31±0.02 


0.48±0.01 


0.74±0.02 


0.56±0.01 


52266.6686" 






2.75±0.04 


0.67±0.02 


1.03±0.02 


0.76±0.02 


52269.6331 






1.42±0.02 


0.44±0.01 


0.74±0.02 


0.68±0.02 


52270.6320 






1.43±0.02 


0.52±0.01 


0.72±0.02 


0.60±0.03 


52271.6485 






1.32±0.03 


0.43±0.01 


0.72±0.02 


0.57±0.01 


52272.6263 






1.15±0.03 


0.44±0.02 




0.42±0.02 


52273.6263 






1.28±0.04 


0.42±0.01 


0.73±0.02 


0.54±0.01 


NOT-SOFIN 


2002 












52508.6869 


2.79±0.14 


0.49±0.01 


1.31±0.02 








52509.6923=" 


3.01±0.18 


0.68±0.02 


1.71±0.03 








52510.7007 


2.48±0.15 


0.49±0.01 


1.29±0.01 








52511.5869 


2.08±0.20 


0.54±0.01 


1.25±0.02 








52512.7255 


2.21±0.20 


0.46±0.01 


1.17±0.01 








52512.7377'' 


1.57±0.18 


0.44±0.01 


1.21±0.02 








52513.7195 


2.06±0.19 


0.44±0.01 


1.29±0.01 








52514.6870 


2.40±0.14 


0.66±0.01 


1.42±0.02 








52515.6157 


2.24±0.18 


0.62±0.02 


1.35±0.02 









^ Flare detected. 

^ Low S/N obtained. 



lines, which are formed at different atmospheric heights. 
Using the spectral subtraction technique described by 
Montes et al. (1995; 1997; 1998; 2000) it is possible to 
study in detail the chromosphere, discriminating between 
different structures: plages, prominences, flares and mi- 
croflares. The synthesized spectra were constructed using 
the program STARMOD developed at Penn State University 
(Barden 1985) and modified by us. The inactive stars used 
as reference stars in the spectral subtraction were observed 
during the same observing run than the active stars (listed 
in Table nj. Representative spectra in the Ca 11 H & K, 
H5, H7, H/3, Ha and Ca 11 IRT (A8498, A8542) line re- 
gions of the quiescent state of PW And are presented in 



Fig. ^ For each region we have plotted the observed spec- 
trum (solid-line) and the synthesized spectrum (dashed- 
line) in the left panel and the subtracted spectrum in the 
right panel. Ha emission above the continuum is detected 
showing an autoabsorption feature. Filled-in absorption in 
other Balmer lines is also detected, and strong emission is 
observed in Ca 11 H & K, H<: and Ca 11 IRT lines. 



The equivalent width {EW) (measured in the sub- 
tracted spectra) for the Ca 11 H & K, He, US, H7, H/3, 
Ha, and Ca 11 IRT (A8498, A8542, A8662) lines for each 
observation is given in Table 01 An estimate of the error 
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Table 4. Absolute surface flux of the different chromospheric activity indicators 
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6.544 
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2.65 


3.33 


1.32 


51387.0396 
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6.043 
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6.111 
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3.52 
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6.441 
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5.561 
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6.493 


6.043 


6.156 
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3.15 
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2.71 




52515.6157 
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" Flare detected. 

^ Low S/N obtained. 



for the EW can be obtained from the relation given by 
Cayrel (1988): 



AEWx ^ 



S/N 



(1) 



where uj is the FWHM of the line, Sx the pixel size in A, 
and S/N the signal-to-noise ratio per pixel in the contin- 
uum. With a typical uj of 1.35 A for the Hq line, Sx — 0.09 
A and S/N ~ 100 for our observations, the error estimate 
is in the order of 6 mA. This equation does not take into 
account errors in the normalization, as well as possible 
blends, and must be used with the observed spectrum in- 
stead of the subtracted one. We have estimated the errors 



in the measured EW (see Table|3l taking into account the 
rms error obtained from the fit between observed and syn- 
thesized spectra in the regions outside the chromospheric 
features and the standard deviations resulting from the 
EW measurements. Deviations in the normalization, dif- 
ferences between the active and the reference stars, and 
internal precisions of STARMOD (0.5 - 2 km s~^ in velocity 
shifts, and ± 5 km in v sini) are taken into account by 
this method. Errors in the chromospheric features in the 
blue spectral region are larger due to the lower S /N of the 
spectra in this region. As an indication of the accuracy of 
the data, we give in Tabled the S/N in the Ca ii H & K, 
and Ha line regions. 
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Ha He I & Na I D1,D2 



Night 7 



Subtracted 




6560 6565 
ITavelenglh (A) 



Wavelength (A) 



Fig. 5. Comparison between pre-flare (night 3), flare max- 
imum (night 4) and gradual decay (night 7) in Ha line in 
the observing run HET-HRS 2001/12. Dashed lines on the 
left panel correspond to the reference star. Dashed and 
dotted-dashed lines on the right panel are for broad (B) 
and narrow (N) components respectively. 




5660 5890 5S00 

Wavelength (A) 



5BB0 5690 5900 

Wavelength (A) 



Fig. 6. Comparison between pre-flare (night 3), flare max- 
imum (night 4) and gradual decay (night 7) in He i D3 
and Na i (Di,D2) lines which are present as emission in 
the subtracted spectrum at flare maximum in the observ- 
ing run HET-HRS 2001/12. Dashed hues on the left panel 
correspond to the reference star. 



The EW have been converted to absolute chromo- 
spheric flux at the stellar surface by using the calibration 
of HaU (1996) as a function of {B-V). In TableHwe give 
the absolute flux at the stellar surface (logics) for the lines 
listed in Table El 

In Table 21 we also give the ratio of excess emission 
EW, ^^(H^j ; for all the observations of PW And, and 
the ratio of the excess emission .^smi a,nd with the 



correction: 



(2) 



given by Hall & Ramsey (1992) that takes into account 
the absolute flux density in these lines and the color dif- 
ference in the components. 

We have used this ratio as a diagnostic for discrimi- 
nating between the presence of plages and prominences 
on the stellar surface, following the results of Hall & 
Ramsey (1992) who found that low -Ehq/^h^ {"^ 1-2) 
can be achieved both in plages and prominences viewed 
against the disk, but that high ratios 3-15) can only 
be achieved in extended regions viewed off the limb. The 
high ratio {Ena/ En^yS, see Table O that we have found 
in PW And indicates that the emission would arise from 
extended regions (prominences). £'8542 /£'8498 ratios arc in 
the range ^ 1-2 which is indicative of optically thick emis- 
sion in plage- like regions, in contrast with the prominence- 
like material inferred by the -Ehq/^'H/3 ratios. 



5. Flare state 

Optical flares are commonly observed in dMe stars, how- 
ever, in more luminous stars flares are usually only de- 
tected through UV or X-ray observations (e.g., Landini 
et al. 1986; Hiinsch & Rcimers 1995; Ayres et al. 1994). 
Optical flares are rare for K-type stars and only a few of 
them have been reported. Monies et al. (1999) detected 
a strong optical flare on LQ Hya, a K2 dwarf star very 
similar to PW And, with a duration of at least 5 hours. 

Abbott et al. (1995) reported a photometric flare on 
PW And in the U band (AC/ ~ 0.2). Two optical flares 
are detected on PW And in our data set. The first one 
was observed during the fourth night of the HET/HRS 
12/2001 observing run (2001 December 23), and exhib- 
ited an enhancement in both Ha and the Ca 11 IRT lines. 
The second flare was observed during the NOT-SOFIN 
2002/08 observing run (2002 August 23) and while it ap- 
pears to be a less powerful event it is possible that the 
flare maximum may have occurred between observations. 

5.1. Description of the flare on 2001 December 23 

The H I Balmer lines are in emission in the observed and 
subtracted spectra at flare maximum and quiescent state 
as can be readily seen in Fig. O . On the other hand, 
Na I Di and D2 lines are clearly revealed only in the sub- 
tracted spectra. The usual chromospheric flare diagnostic 
lines He i D3 and the Mg i b triplet are only observed at 
flare maximum. Observed and subtracted spectra in the 
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Fig. 7. Comparison between pre-flare (night 1), flare max- 
imum (night 2) and quiescent (night 3) in Ha line in the 
observing run NOT-SOFIN 2002/08. Dashed lines on the 
left panel correspond to the reference star. Dashed and 
dotted-dashed lines on the right panel are for broad (B) 
and narrow (N) components respectively. 



Night 3 
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Wavelength (A) 
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Fig. 8. Comparison between pre-flare (night 1), flare max- 
imum (night 2) and quiescent (night 3) in He i D3 and 
Na I (Di,D2) lines which are present as emission in the 
subtracted spectrum at flare maximum in the observing 
run NOT-SOFIN 2002/08. Dashed lines on the left panel 
correspond to the reference star. 



region of He i D3 and Na i doublet are plotted in Fig.Elfor 
the quiescent state (night 3, 2001 Dec 22), flare maximum 
(night 4, 2001 Dec 23) and gradual decay (night 7, 2001 
Dec 26). 

The Ha emission equivalent width (EW) in the sub- 
tracted spectra, increases a factor of ~ 2.1 from the third 
to the fourth nights of the observing run (see Fig. and 
Table EJ, while the Ca 11 IRT lines increases a factor of 
-1.4 (see Table 13 . In Table we list the measured EW 
of He I D3, Na i doublet and Mg i b triplet in the sub- 
tracted spectra for the flare spectra and the presumed qui- 
escent spectra. Emission in the Na i doublet is observed in 
the subtracted spectra throughout the flare event, being 
stronger at flare maximum, (see Fig. while emission in 
He I D3 and Mg i b triplet is present only at the observed 
flare maximum. The detection of prominent He i D3 emis- 
sion indicates that we are observing an energetic flare in 
PW And, as this feature is typically only observed in flare 
events in very active stars such as RS CVn systems and 
UV Ceti flare stars (see Monies et al. 1999 and references 
therein). The Fe 11 A5169 A hne near the Mg 11 b triplet 
is also observed in emission at the flare maximum (see 
Table El . 

One notable feature is the existence of a broad emission 
component in both the subtracted and observed profiles in 
Ha at the flare maximum. A two Gaussian fit is necessary 
to obtain a good match to the total emission: a narrow 
component (N) having a FWHM of — 59 km s^^ and a 
broad component (B) with FWHM ~ 242 km s~^ The 
profile is also asymmetric with the B component blue- 



shifted. This could be interpreted as high velocity mass 
ejection (Houdebine et al. 1990) or high velocity chromo- 
spheric evaporation (Gunn et al. 1994). A broad emission 
component is also detected in the Ha subtracted profile 
during the gradual decay, yielding a FWHM of 53 km 
and 115 km s^^ respectively for the N and B components. 
Again, the profile is asymmetric with the B component 
showing a slight blue-shift. The presence of a red-shifted B 
component in the subtracted profile in night 3 suggest that 
the flare event might have begun this night, taking — 24 
hours to reach maximum (night 4) . This red asymmetry is 
often interpreted as the result of chromospheric downward 
condensations (CDC) (Canfield et al. 1990 and references 
therein) . No significant B component is observed in quies- 
cent state. The EW{)L^.a) returns on night 10 to the same 
value as observed in the quiescent state in night 2 (see 
Table EJ . Analyzing this EW data, we estimate a dura- 
tion of approximately six days for this flare event. 

In TableElwe fist the intensity (I), FWHM and EW of 
the B and N components and the shift between them (An - 
Ab), as well as the contribution of B and N components 
to the total EW in the Ha line. The total values of I, EW 
and logF are also given. 

5.2. Description of the flare on 2002 August 23 

The second flare occurred during the NOT-SOFIN 
2002/08 observing run and appears to be a less energetic 
or more rapid event. The excess Ha emission EW in the 
subtracted spectra increases in a factor of ~ 1.46 from the 
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Table 5. Measured EW{A) of He i D3, Na i Di and D2 and the Mg i b triplet in the subtracted spectra during the 
flares. 



Date 


UT 


He I D3 


Na I Di 


Na I D2 


Mg I bi 


Mg I b2 


Mg I bs 


Fo II A5169 A 


HET-HRS 2001/12 


22 Dec (pre-fiare) 

23 Dec (flare max) 

26 Dec (gradual decay) 


03:44 
04:00 
03:09 


0.208±0.006 


0.087±0.010 
0.255±0.006 
0.074±0.010 


0.123±0.012 
0.273±0.008 
0.102±0.014 


0.042±0.001 


0.045±0.001 


0.062±0.001 


0.082±0.001 


NOT-SOFIN 2002/08 


22 Aug (pre-flare) 

23 Aug (flare max) 

24 Aug (quiescent) 


04:16 
04:24 
04:39 


0.023±0.001 


0.099±0.002 
0.112±0.002 
0.093±0.001 


0.147±0.002 
0.167±0.002 
0.147±0.001 











Table 6. Ha line parameters (total and narrow and broad components) measured in the subtracted spectra during 
the flares 







Hq Total emission (T) 


Hq 


Broad component (B) 


Ha Narrow com 


iponent (N) 




Date 


UT 


I 


EW 


logF 


I 


FWHM 


EW 


B/T 


I 


FWHM 


EW 


N/T 


AA 








(A) 






(A) 


(A) 


% 




(A) 


(A) 


% 


An - -^B 


HET-HRS 2001/12 


21 Dec (quiescent) 


03:37 


0.75 


1.16 


6.478 




















22 Dec (pre-flare) 


03:44 


0.82 


1.31 


6.531 


0.13 


2.94 


0.41 


31 


0.72 


1.18 


0.90 


69 


-0.14 


23 Dec (flare max) 


04:00 


1.06 


2.75 


6.853 


0.29 


5.29 


1.66 


60 


0.79 


1.30 


1.09 


40 


0.10 


26 Dec (gradual decay) 


03:09 


0.81 


1.42 


6.566 


0.14 


2.52 


0.37 


26 


0.69 


1.16 


1.05 


74 


0.04 


NOT-SOFIN 2002/08 


22 Aug (pre-flare) 


04:16 


0.83 


1.31 


6.531 


0.06 


4.98 


0.31 


24 


0.80 


1.23 


1.00 


76 


0.76 


23 Aug (flare max) 


04:24 


0.96 


1.71 


6.647 


0.09 


5.16 


0.48 


28 


0.90 


1.28 


1.23 


72 


-0.40 


24 Aug (quiescent) 


04:39 


0.81 


1.29 


6.524 





















quiescent state to the flare maximum (significantly differ- 
ent from the 2.1 factor obtained in HET-HRS 2001/12). 

In Table El we list the measured EW of He i D3 and 
Na I doublet in the subtracted spectra. The Mg i b triplet 
is not observed at this epoch due to a wavelength gap in 
the configuration of the spectrograph. Emission in the Na i 
doublet is observed in the subtracted spectra in both flare 
spectra and is stronger at flare maximum. The emission in 
He I D3 is present only in the flare maximum spectrum (see 
Fig.|SI). The much lower EW{lle i D3) and the increment 
in ^^^^(Na i DI & D2) measured at flare maximum both 
indicate a less energetic event than that observed on 2001 
December 23. 

As in the case of HET-HRS 2001/12 run, a broad emis- 
sion component (B) in both the subtracted and observed 
Ha profile is present at flare maximum (see Fig. d) . The 
two-Gaussian fit used yields a FWHM ~ 236 km s~^ for 
the B component and FWHM ^ 58 km s^^ for the nar- 
row one (N). In the pre-flare phase we obtain a FWHM 
for the N component of 56 km s~^, very similar to that 
obtained for N component in HET-HRS 2001/12 (53 km 
s-i). However, the FWHM of 228 km s^^ in the B com- 
ponent is very different of 134 km s^^ from HET-HRS 
2001/12. No significant B component is detected on the 
night 3; 24 Aug. This suggests that flare maximum likely 
occurred between night 1 and night 2 and that the max- 



imum observed on night 2 is part of the gradual decay 
phase. Thus, the duration of this event is between 2 and 
3 days, shorter than in 2001/12. Contrary to what is ob- 
served in the HET-HRS 2001/12 event, the B component 
is blue-shifted in the pre-flare observation and red-shifted 
at the maximum (night 2). Different dynamic processes 
could be taking place during this Hare event. Another pos- 
sible explanation is that these asymmetries could be due 
the position of the active region over the disk. 

In Table El we list the Ha hue parameters I, FWHM, 
EW, and relative contribution of the B and N compo- 
nents, as well as the shift between them (An - Ab) during 
the flare. The total values of I, EW, and logF of the total 
excess Ha emission are also presented. 



6. Photospheric variations 

Since PW And is a young K2-dwarf with a rotational ve- 
locity of 22.6 km s^"'^, it is expected to show large dark 
photospheric spots which produce rotational modulation 
of the light coming from its surface. A photometric period 
was obtained by Hooten & Hall (1990) by fitting a sine 
curve to the data using a standard less squares technique. 
Values of 1.64 ± 0.01, 29 ± 1 and 1.745 ± 0.005 days were 
calculated for three different epochs. They propose 1.745 
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Table 7. Heliocentric radial velocities and CCF bisector 
parameters 



Fig. 9. Left: Example of a cross-correlation function 
(CCF). Abscissa is given in units of pixels. The peak of 
the CCF has been fitted by a cubic-spline. Right: A zoom 
of the bisector is plotted in the right panel. 

days as the more reliable period is it had the least internal 
scatter. 

To study spectroscopic variations produced by spots 
we have used the fact that cool spots moving across the 
disk of a star will produce a weaker contribution to the 
total integrated line profile. Thus the spectra -integrated 
over the disk will show a bump traveling across the pro- 
file as the star rotates. Two methods are commonly used 
to interpret this information. Doppler imaging was devel- 
oped with the aim of obtaining a two dimensional picture 
of spots (Vogt & Penrod 1983, Vogt et al. 1987, Collier- 
Cameron & Unruh 1994, Rice et al. 1989, Strassmeier & 
Rice 1989). Here high signal-to-noise ratio (S/N) and res- 
olution (A/ A A ~ 100000) are needed in order to obtain 
both accurate position and dimensional information on 
the spot(s). Another method, first used in the case of the 
Sun, can be adapted to stars using high A/ A A and S/N 
spectra (Toner & Gray 1988). The technique consists of 
calculating a bisector for line profile, i.e. the middle points 
of the line profile taking points of equal intensity in both 
sides of the line. Variations of the bisector are related to 
the existence of surface features moving across the disk, as 
seen in ^ Boo A by Toner & Gray (1988) and in A And by 
Donati et al. (1995). Both techniques have serious limita- 
tions (Dempsey 1991). Doppler imaging requires the star 
to be rotating rapidly wsini > 36 km s~^ to give suffi- 
cient spatial resolution. Values of i over 10-30 degrees are 
needed to guarantee high rotational broadening. However, 
the distortions induced by spots are of the order of 1-10 
per cent of the continuum, and S /N > 300 is required to 
detect them. Moreover, bisector analysis is only appropri- 
ate for very sharp lines, vsini < 10 km s~^, restricting 
the candidates to the very slow rotators or stars with very 
low value of inclination. S/N > 150 is also needed to apply 
this technique. A common problem with both techniques 
is the use of only one unblended line for the study, al- 
though several lines can be combined in order to obtain 
more accurate results. 

A powerful method was developed and applied to sev- 
eral stars in the range of 15 < wsini < 40 km s~^ 
by Dempsey et al. (1992) using a correlative analysis. A 



Date 


UT 


Vhol 






Phasc*^ 




(hh:mm) 


(km s"i) 


(A) 


(km s"i) 




2.2m-FOCES 1999 


25 Jul 


04:09 


-11.76 ± 0.59 


0.043 


1.90 


0.00 


26 " 


00:57 


-11.49 ± 0.23 


-0.007 


-0.35 


0.49 


27 " 


02:25 


-12.50 ± 0.37 


0.070 


3.14 


0.10 


28 " 


00:57 


-10.51 ± 0.42 


-0.044 


-1.97 


0.64 


29 " 


03:01 


-11.95 ± 0.26 


0.040 


1.78 


0.26 


30 " 


01:57 


-11.11 ± 0.47 


0.037 


1.63 


0.81 


HET-HRS 2001-2002 


20 Dec 


04:12 


-13.08 ± 0.41 


0.038 


1.72 


0.73 


21 " 


03:37 


-13.77 ± 0.30 


0.014 


0.65 


0.29 


22 " 


03:44 


-12.45 ± 0.35 


-0.022 


-0.99 


0.86 


23 " 


04:00 


-11.51 ± 0.34 


-0.043 


-1.96 


0.44 


26 " 


03:09 


-12.40 ± 0.40 


0.046 


2.07 


0.14 


27 " 


03:07 


-14.83 ± 0.35 


0.043 


1.95 


0.71 


28 " 


03:31 


-12.34 ± 0.30 


-0.014 


-0.64 


0.29 


29 " 


02:52 


-15.04 ± 0.46 


-0.001 


-0.07 


0.85 


30 " 


02:59 


-11.90 ± 0.33 


-0.038 


-1.71 


0.43 


NOT-SOFIN 2002 


22 Aug 


04:16 


-11.62 ± 0.28 


0.047 


2.09 


0.13 


23 " 


04:24 


-11.65 ± 0.34 


0.042 


1.86 


0.71 


24 " 


04:39 


-11.73 ± 0.25 


0.055 


2.44 


0.29 


25 " 


01:55 


-11.96 ± 0.33 


0.108 


4.84 


0.79 


26 " 


05:15 


-8.83 ± 0.26 


-0.054 


-2.43 


0.45 


26 " 


05:37 


-8.81 ± 0.32 


-0.034 


-2.15 


0.45 


27 " 


05:06 


-9.95 ± 0.29 


0.014 


0.60 


0.02 


28 " 


04:19 


-9.03 ± 0.40 


-0.031 


-1.42 


0.57 


29 " 


02:37 


-10.61 ± 0.34 


0.015 


0.67 


0.10 



" fphot = 1.745 days (from Hooten & Hall 1990). 

non active star (template) is cross-correlated with the ac- 
tive star producing a cross-correlation function (CCF). 
Variations in the peak of the CCF are related to changes in 
the line profiles caused by spots. To quantify the temporal 
variations in the CCF the bisector of the peak of the CCF 
can be calculated. This CCF bisector is not the same as 
the bisector described by Toner & Gray (1988), since the 
bisector of a line measures velocity fields while the CCF 
bisector quantifies the asymmetry of the CCF. The best 
results are obtained when both template and active stars 
have a similar spectral type. The advantage of this tech- 
nique is the possibility of using many absorption lines for 
the calculation of the CCF (Dempsey et al. 1992). As the 
information of the lines is redundant, there is less restric- 
tion in the S/N (-150). Resolution of A/AA - 40000 is 
sufficient to obtain accurate results. We employ this latter 
technique to study the effects of star-spots on the profiles 
of photosphcric lines in PW And. 

Three observing runs with the most extensive rota- 
tional phase coverage are used to study the effects of 
photosphcric spots. These are the 2.2m-FOCES 1999/07, 
HET-HRS 2001/12 and NOT-SOFIN 2002/08 observa- 
tions (see Section 2 and Tabic ^J. The CCFs were deter- 
mined using the IRAF routine FXCOR in the same proce- 
dure as used in the determination of the radial velocities 
Vhol (see Section 3.2). Here we limit calculations to the 
wavelength regions ranging from 6300 to 6465 A and 6670 
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to 6760 A (see Fig.lllland I16|l . which includes lines com- 
monly used in Doppler imaging like Fe I lines in 6411.644 
and 6430.841 A and Ca i lines in 6439.073, 6462.566 and 
6717.680 A. The asymmetries of the CCF are very subtle 
in PW And necessitating calculation of its bisector. The 
determination of the CCF bisector has been carried out 
in the usual manner, taking points of equal intensity for 
the calculation of the middle points. For each point on the 
left of the CCFs, a matching point is found on the right 
using a cubic-spline interpolation (see Fig. EJ. In order 
to quantify the changes in the CCF bisector, the differ- 
ence between the mean of ten interpolated points in the 
top and in the base of the CCF bisector (AAbis) has been 
measured. We have choose a fixed value of 0.2 as the base 
of the CCF (see Fig.^ to avoid the noise at the zero level 
of the CCF. In the subsequent discussion the values of 
AAbis have been converted into velocities (Awb) for better 
understanding the science. 

We have analyzed the CCF of each observing run sepa- 
rately, in order to avoid the possible slightly differences of 
the instrumental profile of the different spectrographs, and 
because the behaviour of the photospheric activity could 
be different at these three epochs separated by nearly 
three years. 

6.1. CCF bisectors in tlie observing run 2.2m-F0CES 
1999/07 

The six observations over six nights cover almost three 
rotational cycles. GJ 706 (K2 V) has been used as tem- 
plate to calculate the CCF for PW And. In Fig. Uni we 
have plotted in phase the CCF bisectors and the values of 
Awb using a period of P = 1.745 days from Hooten & Hall 
(1990) with the time of the first observation (To) taken as 
phase 0.0. 

The maximum Aub measured is 3.14 km s~^ (see 
Table 13) at phase 0.1, changing to -1.97 km at phase 
0.64. The sequence is compatible with the existence of two 
spots at nearly opposite longitudes, one of them near the 
pole since a perturbation in the core of the line is present 
in every phase. In Fig. ^2 we have plotted the profile of 
Fe I 6430.84 A and Ca i 6439.07 and 6717.68 A photo- 
spheric lines. A disturbance in the core of the Ca I lines, 
slightly shifted with phase, is clearly present for all phases. 

The variations observed in the radial velocity (T4iei), 
determined as is explained in Section 3.2, are related to the 
existence of spots on the star. A clear correlation between 
Avh and Vhei is evident (see Fig. [T^. 

6.2. CCF bisectors in tiie observing run HET-HRS 
2001/12 

Nine observations of PW And were obtained during this 
observing run covering 6.3 rotational periods with a gap of 
two days between nights 4 and 7. In order to obtain accu- 
rate results for CCF bisectors, the non-active star GJ 706 
from 2.2m-FOCES 1999/07 was used as template. The 
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Fig. 10. CCF bisectors (lower panel) and Aub (upper 
panel) for the observations of 2.2m-FOCES 1999/07 run 
arranged in phase using the photometric period given by 
Hooten & Hall (1990) and the first observation as phase 
0.0. 



Fe I (X6430.84) Ca I (X6439.07) Ca 1 (X6717.68) 




6429 6430 6431 6437 6438 6439 6716 6717 6718 



Wavelength (A) 

Fig. 11. Line profiles (AA6430.84, 6439.07 and 6717.68 A) 
of PW And obtained during 2.2m-FOCES 1999/07. Phase 
is given at the left panel, for each observation. The same 
pattern is seen in the three lines for each phase. 
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• rOCES 1999/07 

• HET 2001/12 
A NOT 2002/08 




-8 -10 -12 -14 -16 

Radial velocity (km s ^) 

Fig. 12. Avh in the bisector vs. Vhci- Changes in T4oi are 
related to the existence of spots on the star. Differences in 
T4ioi between the different observing runs could be due to 
the different radial velocity standard stars used. Changes 
in the slope of the lines could be due to changes in the 
pattern of spots. 



spectrum of GJ 706 was convolved with a Gaussian in or- 
der to degrade the spectral resolution to the lower A/ A A 
in HET-HRS 2001/12. 

CCF bisectors have been plotted in Fig. ^1 using the 
same To and Pphot as in 2.2m-FOCES 1999/07 with the 
aim of comparing both observing runs. In spite of the more 
than two years time difference, the change in the CCF 
bisector follows a sequence very similar to that for 2.2m- 
FOCES 1999/07 when arranged in phase. This could be 
due to the existence of large patterns of spots persisting 
with time, perhaps due to active longitudes. 

As in the previous observing run, a correlation between 
Awb and Viici is found (see Fig ll2l) . In this case the slope is 
different and, due to the lower spectral resolution of this 
run, a larger scatter in Vhci is observed. 

6.3. CCF bisectors in tlie observing run NOT-SOFIN 
2002/08 

For this data set we used the template K2-dwarf star HR 
222 (observed during every night) as both the radial ve- 
locity and spectral type standard in calculating the CCF. 
Due to a large gap in the spectral region between 6420 and 
6540 A caused by the configuration used here, a range 
from 6360 to 6415 A has been used for measuring CCF 
bisectors instead of that employed in the previous observ- 
ing runs. In this region the strong unblended Fe i lines 
AA 6393.72, 6408.03 and 6411.64 A are present. The Fe i 
lines AA 6400.33 and 6400.01 are blended and have been 
rejected for this analysis (see Fig. I14|l. 

As in the preceding cases, the photometric period of 
1.745 days has been used to calculate the phase of the 




Phase 

Fig. 13. CCF bisectors and Aub for the observations of 
HET-HRS 2001/12 arranged in phase using the photo- 
metric period given by Hooten & Hall (1990) and To of 
2.2m-FOCES 1999/07 as phase 0.0. 





Fel 6393.72 


Fel 6A 


08.03 Fel 6411.64 - 


■ PW 


And ' 






: HR 


222 








6390 


6395 


6400 6405 
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6410 



Fig. 14. Observed spectrum of PW And and HR 222 in 
the region between 6360 to 6415 A used in the determi- 
nation of the CCF in NOT-SOFIN 2002/08. The blended 
AA 6400.33 and 6400.01 have been rejected in this work. 



CCF bisectors. The results have been plotted using Tq of 
2.2m-FOCES 1999/07 as phase 0.0 (Fig.U^). 

A clear difference in the phase between this and the 
previous observing runs is present. Nevertheless, the se- 
quence is similar to that in HET-HRS 2001/12 and 2.2m- 
FOCES 1999/07. The maximum Avb measured is 4.84 km 
s"""^ (see TableUJ at phase 0.79. This is larger than the ob- 
tained on July 1999, and less than 2.50 km s^"'^ is measured 
for the rest of the phases. This suggests the presence of a 
large spot on the approaching limb. 

The spot disturbance in the line profile is not the same 
as in 2.2m-FOCES 1999/07. The presence of a perturba- 
tion in the core of absorption lines is not clear in Fig. 1161 
although a feature disturbing the profile is shown. There 
is no apparent evidence for a polar spot. 

The correlation between Awb and Vhci is plotted in 
FigEl 111 this case the slope of the line is similar to that 
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Fe I (X6393.72) Fe 1 (X6408.03) 



o 
o 



W 



Phase 



Fig. 15. CCF bisectors and Avh for the observations of 
NOT-SOFIN 2002/08 arranged in phase using the photo- 
metric period given by Hooten & Hall (1990) and To of 
2.2m-FOCES 1999/07 as phase 0.0. 



of HET-HRS 2001/12 observing run. These changes in the 
slope from one observing run to another could indicate 
changes in the pattern of the spots at these epochs. 

7. Relation between photospheric and 
chromospheric activity 

When studying the chromospheric activity indicators in 
Section 4, variations in the EW of several lines, such as 
Ha, H/3 and Ca ii H & K were found (see Tables 01 & 13 . 
When arranged in phase, the EW(Ha) and EW{Ca ii H 
& K) appear to follow a clear sequence for each observ- 
ing run. Comparison with results from the CCF bisectors 
shows a correlation between photospheric and chromo- 
spheric activity. In Figs. El El & El we have over- plotted 
AAbis and Vhei with the EW of chromospheric lines for 
each observing run. 



Fe 1 (X6411.64) 




Wavelength (A) 

Fig. 16. Line profiles (AA 6393.72, 6408.03 and 6411.64 A) 
of PW And obtained during NOT-SOFIN 2002/08. Phase 
is given at the left panel, for each observation. 



jecture is further supported by the values of £"3542 /^'8498 
calculated in Section 5 and listed in Table 01 



7.1. 2.2m-F0CES 1999/07 

For this run we study the chromospheric lines Ha and 
Ca II H & K in relation to the photospheric line behavior 
(see Fig.ini- The EW(G'a ii H) and EW{G& 11 K) have 
been summed together due to the low S/N of the contin- 
uum in this region of our spectra (see Table Note that 
the EW of Ca 11 lines could not be measured on night 4, 
which corresponds to phase 0.64. 

A clear relation between £'W^(Ca 11 H + K) and the 
variations in Awb and Vjici is observed (see Fig. ll7|l . For the 
Ha line, the relation is not as clear, but a similar trend is 
followed from phase 0.5 to 1.0. This supports the existence 
of active regions connected to the photospheric features. 
The presence of plage-like regions can be inferred from 
the variation on Ca 11 H & K lines emission and this con- 



7.2. HET-HRS 2001/12 

The chromospheric lines Ha and Ca 11 A8542 A are used 
for the comparison during this run since Ca 11 H & K lines 
were beyond wavelength range observed in this epoch (see 
Section 2). 

In Fig.Elwe have plotted the EW of these lines as well 
as CCF bisectors and Vhci- The Ca 11 A8542 A line does not 
appear to be a good indicator of changes in the activity of 
the star at this epoch. Only the flare on 2001 December 
23, corresponding to phase 0.44, is clearly visible using this 
line. Taking into account that Vhci is a minimum at this 
phase, the flare occurred when the photospheric feature 
was near the limb. 

Variations in Ha follow a pattern similar to that found 
in CCF bisectors and Vhoi only up to phase 0.5. The less 
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Fig. 17. Comparison between photospheric (Awb, Vhci) 
and chromospheric {EW(Ha), EW {Ca. ii H & K)) varia- 
tions in 2.2m-FOCES 1999/07. 



more patchy phase sampling in this data set makes any 
correlation less clear. 



7.3. NOT-SOFIN 2002/08 

For this observing run we have again the Ha and Ca ii K 
lines. The flare on 2003 August 23 corresponds to phase 
0.71 and is clearly visible in both chromospheric lines (see 
Fig.Uni). Contrary to the event in HET-HRS 2001/12, the 
flare does not appear to have occurred over the limb, but 
seems more associated with an active region on the disk. 
An overall correlation between both Ha and Ca ii K lines 
with CCF bisectors and Vhc\ is apparent. This suggests 
that the chromospheric active regions are connected to 
the photospheric features in the same way as 2.2m-FOCES 
1999/07, although the correlation of Ha is perhaps more 
clear at this epoch. 

In order to look for any features in the Ha emission 
line that are related to the variations seen in the photo- 
spheric lines, we have applied the bisectors method to the 
emission profile in the subtracted spectrum (see Fig. I21|) . 
Again, for each point on the left of the Ha subtracted pro- 
file, a matching point on the right is found with a cubic- 
spline interpolation to obtain a bisector of the subtracted 
emission line. Any features moving across the profile must 
disturb this bisector in the same way as in the case of 
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Fig. 18. Comparison between photospheric and chromo- 
spheric variations in HET-HRS 2001/12. 



a CCF bisector. The resulting profile bisectors and Avh 
are plotted in Fig. 1201 using the same Pphot and Tq as in 
Fig. The scale of the bisectors has been changed to 
make it easier to see the curvature. Only the top of the 
profile is perturbed by the features, while the bottom is 
less sensitive, probably due to the difference in flux and 
S/N between the core and the wings of the emission. In 
order to avoid changes in the bottom of the profile bisector 
induced by changes in the broad component, as is present 
in the fiare state, we have discarded ^--^20 per cent of the 
bisector nearest the continuum. The sequence obtained in 
Fig. |2D1 is similar to that of Fig. ^1 suggesting the pres- 
ence of cool chromospheric material associated with the 
photospheric features. This is consistent with the results 
of the EW of chromospheric emission lines (see Fig. I19|l . 



8. Conclusions 

In this paper we have analyzed high resolution echelle ob- 
servations, taken during the period 1999-2002, to describe 
in detail the activity on the young K2-dwarf star PW And. 
We also have presented new determinations of key physi- 
cal parameters. In addition, a study of the lines profile has 
been carried out in order to relate the variations found in 
the heliocentric radial velocity (Viici) and chromospheric 
emission with the existence of photospheric features co- 
rotating with the star. 
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Fig. 19. Comparison between photospheric and chromo- 
spheric variations in NOT-SOFIN 2002/08. 
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Fig. 20. Bisectors of the subtracted emission line profiles 
of Ha for the observations of NOT-SOFIN 2002/08. 



Radial velocity standard stars have been used to calcu- 
late accurate heliocentric radial velocities for each obser- 
vation by using the cross-correlation technique. Variations 
in the Vhei up to 3 km in the same observing run 
were found. The mean value of the velocities of all the 
observing runs (-11.15 ± 0.05 kms~^) is used to calcu- 




6560 6S61 6562 6563 6564 

Wavelength (A) 

Fig. 21. Subtracted Ha spectra for the NOT-SOFIN 
2003/08 observing run ordered in phase. Flare occurred 
for phase 0.71. 



late the Galactic space- velocity components {Uy,W) and 
confirms its membership to the Local Association moving 
group. By using the information provided by the width 
(FWHM) of the peak of the cross-correlation function 
(CCF), we have determined a projected rotational ve- 
locity, wsini, of 22.6 kms~^. The radius for a K2V star 
(i? = 0.80 Rq) is compatible with the minimum radius 
{R = 0.78 ± 0.11 R0) obtained from fsini and the pho- 
tometric period. The age of PW And has been bounded 
between 30 to 80 Myrs from its position in relation to 
pre-main sequence isochrones in the color magnitude di- 
agram. The measured EW{Li I) of 273 mA indicates an 
age similar to that of the stars in the Pleiades cluster or 
younger, which is in agreement with the estimate from the 
color-magnitude diagram. 

Using the spectral subtraction technique, all the opti- 
cal activity indicators from the Ca li H & K to Ca 11 IRT 
lines have been analyzed. Strong emission in the observed 
spectra is present in Ha, Ca 11 H & K and IRT, the re- 
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mainder of the Balmer lines are clearly visible in emission 
in the subtracted spectra. The ratio that we have 

found in 2.2m-FOCES 1999/07 and NOT-SOFIN 2002/08 
suggest that the emission of these lines could arise from 
prominence-like material, whereas the ratio of the excess 
emission -fsMa. indicates that the Ca II IRT emission arise 

-C/8498 

from plagc-like regions. Two flares were detected in differ- 
ent epochs (2001 and 2002) which confirms the presence 
of strong magnetic activity on the star. At the same time, 
variations in EW{B.a) and EW{Ca, II H & K) arc observed 
in all of the observing runs further confirming magnetic 
surface activity. 

A detailed analysis of the profiles of the photospheric 
absorption lines have been done by calculating the bisec- 
tors of the CCF resulting of cross-correlate the spectra of 
PW And with the spectrum of a non active star of similar 
spectral type. These CCF bisectors change with a period 
similar to the photometric period of the star, suggesting 
the presence of features (probably cool dark spots) in the 
photosphere which are co-rotating with the star. The vari- 
ations of the CCF bisectors found in the three epochs an- 
alyzed confirm the prevalence over time of large spots or 
the existence of active longitudes were spots are continu- 
ously being generated. The chromosphcric active regions 
appear to be associated to the photospheric features, as a 
clear correlation between the variations observed in chro- 
mosphcric activity and in the photospheric lines is found. 
Finally, variations in the bisector of the subtracted Ha 
line further support the presence of chromospheric active 
regions, related to the cool photospheric spots. 
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